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An apparatus has been designed and constructed to allow measurements of fugacities in gaseous mixtures 
containing hydrogen. The apparatus makes use of a semipermeable membrane to allow a direct measurement 
of the partial pressure of a permeating component (in this case, hydrogen) in a mixture with a nonpermeating 
component. In this study, measurements were made on mixtures of hydrogen/methane and hydrogen/propane. 

The apparatus is designed to operate at moderate and high temperatures (ambient to 250 °C) and moderate 
and high mixture pressures (3 to 50 MPa). The actual pressure range that is experimentally accessible is 
dependent on the fluid mixture under study, due to its influence on the pressure gradient across the semiper- 
meable membrane. 

The pressure measurements are done isothermally on a series of concentrations of the binary. Each mixture 
is characterized using a developmental gas-chromatograph which has been built especially for this work. The 
gas-chromatograph was calibrated using standard mixtures of hydrogen/methane and hydrogen/propane, pre- 
pared gravimetrically in our laboratory. 

Using measured values of the mixture pressure, hydrogen partial pressure and mixture mole fraction at a 
given temperature, fugacity coefficients were determined using the virial equation. The measured values are 
compared with some previous data and general trends are discussed. 

Key words: fugacity coefficients; fugacity measurements; gas-chromatograph; gas mixtures; hydrogen/methane; 
hydrogen/propane. 



1. Introduction 

The concept of fugacity owes its origin to Lewis 
[1,2]', who desired some measure of "escaping ten- 
dency" upon which to base his idea of equilibrium. His 
original use of partial molar free energy suffered from a 
negative infinity in the case of a gas in the limit of low 
pressure [3]. The concept of fugacity, which was de- 
vised to remedy this difficulty, is a convenient way to 
express the chemical potential of a substance in experi- 
mentally accessible terms. Since defining equations and 
useful functional relationships are presented in a variety 
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of sources [4-6], a discussion of the basics of fugacity 
will not be repeated here. 

Fugacity is a most important quantity in experimental 
thermodynamics because of its theoretical significance 
and general usefulness. The phase equality of tem- 
perature, pressure, and component fugacity is the crite- 
rion for physical equilibrium of a system. In a system 
involving chemical reactions, fugacities appear in the 
equilibrium constants describing the multicomponent 
system. Departures from equilibrium can be described in 
terms of fugacity gradients [7], Since fugacity is a mea- 
sure of nonideality of a gaseous system, measurements of 
fugacity and fugacity coefficients are of value in re- 
search on equations of state, especially work involving 
mixtures. 

There are several techniques by which fugacities in 
gas mixtures may be measured experimentally. The most 
generally used method involves consideration of the 
P-V-T surface of gas mixtures of composition j> s . The 
component fugacities f t are then obtained from [8,9]: 
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\afi=laPXi+f(VyRT-l/P) dP 



(1) 



where V\ is a partial molar volume, and x, is the mole 
fraction of component i, and P, T and R have their usual 
meanings. The P-V-T-x data used in this type of mea- 
surement must be very precise, because of the loss of 
precision of derived thermodynamic functions. In addi- 
tion, large quantities of PVT data are necessary for the 
successful exploitation of the above equation. One can 
also use the above equation with various equations of 
state as a predictive tool for fugacity. 

There are more accurate techniques for experimental 
determination of fugacities in gas mixtures [10]; how- 
ever these are applicable only in specific instances. One 
such case involves mixtures containing hydrogen. For 
these mixtures, a physical equilibrium technique using a 
semipermeable membrane is applicable. The principle of 
this technique is quite simple and is illustrated in figure 
1 . An experimental chamber is divided into two regions 
by a membrane permeable only to component i. The 
fugacity of pure component i at partial pressure Pi will 
be equal to the fugacity of component i at P m , the total 
mixture pressure, after equilibrium has been established. 
The three general criteria for equilibrium are thus satis- 
fied: the equality of temperature, pressure, and fugacity. 
The permeation through the membrane, and hence the 
approach to equilibrium is driven by the equalization of 
chemical potential on both sides of the membrane. 
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Figure 1-Principle of the semipermeable membrane technique. 



The observance and measurement of hydrogen partial 
pressure was originally done using a semi-permeable 
membrane by Ramsay [11] and Tsakalotos [12]. The first 
application of this technique to measurement of gas mix- 
ture component fugacities was due to Krishnamurty [13] 
at Columbia University. His original study of hydrogen- 
helium mixtures was followed by several studies of 
hydrogen-propane mixtures [14, 15], and a few studies of 
ternary mixtures containing butane and ammonia with 
hydrogen-propane [15,18]. The present apparatus differs 
markedly from the original Columbia University instru- 
ment, with many of the difficulties [19,20] encountered 



with the Columbia design being addressed and de- 
scribed in the experimental section. 

Once a hydrogen partial pressure has been measured 
across the semipermeable membrane at a particular tem- 
perature, the fugacity of this single component can be 
determined. It is more convenient to deal with a dimen- 
sionless quantity called the fugacity coefficient, fa, de- 
fined by 



fa=f/(xPd 



(2) 



where P, is the partial pressure of component i [6]. The 
fugacity coefficient of this single component can be de- 
termined using the virial equation (truncated after the 
third virial coefficient) [17,18], 

ln^=f(f) + ^-(^} 0) 

The above quantities are for pure hydrogen only. To 
determine the fugacity coefficient of hydrogen in the 
mixture, we use the defining equation for the fugacity 
coefficient [10]: 



K=/h/W, 



2 * m 



(4) 



where ;c Hl is the mole fraction of hydrogen in the mix- 
ture, and P m is the total pressure of the mixture. The 
value of P m is necessarily greater than that of P it since it 
also includes the partial pressure of the nonpermeating 
component. 

2. Test Systems 

The systems chosen for the initial studies on this in- 
strument are gaseous binary mixtures of hydrogen/ 
methane and hydrogen/propane. These were chosen for 
several reasons. First, there was no evidence that either 
mixture would be detrimental to the membrane, or that 
temporary poisoning of the membrane would cause 
slow diffusion rates. Second, preliminary predictions 
made using the Redlich-Kwong equation indicated that 
these systems would be quite informative [8], 
Hydrogen/propane is expected to have appreciable 
nonideality [15], while hydrogen/methane, on the other 
hand, would be expected to show a smaller change of 
fugacity coefficient with mole fraction. This second sys- 
tem, therefore, provides a test of the overall sensitivity 
of the apparatus, since the total expected change is more 
comparable to the experimental error. Third, there are 
previous data on hydrogen/propane with which to 
compare the present data [18,19]. 

Hydrogen/methane and hydrogen/propane mixtures 
pose no unfavorable analytical requirements. The chro- 
matographic separations are easily done using one 



128 



column for each mixture. Heart-cutting or multi- 
dimensional techniques are unnecessary. The retention 
times are quite reasonable, and well-shaped peaks are 
obtained at experimentally convenient column and in- 
jector temperatures. The analytical conditions for both 
mixtures are thus very conducive to precise quan- 
titation. 

3. Experimental 

3.1 Pressure Vessel 

A schematic diagram of the essential features of the 
apparatus is given in figure 2. The heart of the apparatus 



is a large, thick-walled pressure vessel which serves as 
the experimental chamber. The vessel has a relatively 
large internal volume (2050 ml), so that the removal of 
small aliquots of mixture for analysis will cause minimal 
disturbance of equilibrium. The opening of the vessel is 
also large (11.43 cm diameter), to accommodate the 
membrane manifold, a mixer, and required feed- 
throughs. The vessel was machined from a section of 
work hardened 316 stainless steel barstock, with final 
dimensions and clearances exceeding the requirements 
set forth by the ASME [21,22]. A bolted closure is em- 
ployed in the vessel, using eight bolts machined from 
4340 (AISI designation) steel. The high loading on the 



Vacuum 



H 2 Line tx> 



Exhaust 




Purge 



To 

Secondary 

Sampling 

System 



Oven 



Figure 2-Schematic diagram of the apparatus. 
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11.43 cm diameter closure required the use of an ultra- 
high strength steel for the bolts. Sealing is provided by 
either a 25% glass-filled PTFE gasket or 316 stainless 
steel gasket (the stainless steel gasket is used for tem- 
peratures above 220 °C). The vessel with its closure is 
capable of containing a pressure of 50 MPa at 340 °C. 
The lid of the vessel also accommodates a thermometer 
probe well, and a provision for a safety head (equipped 
with a 41 MPa rupture disk) and supports a propeller- 
type mixer. The mixer is a commercial magnetic linkage 
type with air cooled magnets. 

3.2 Membrane Manifold 

Due to the importance of the membrane manifold for 
this experiment, and its attendant complexities of fabri- 
cation, a detailed discussion is provided. The semiper- 
meable membrane manifold, pictured in figure 3, is sus- 
pended from the underside of the vessel lid using a 
compression Fitting. 

As stated previously the working principle of the 
semi-permeable membrane method depends upon the 
separation of the mixture (i, j, . . .) and a pure, perme- 
ating component (i) into two separate "compartments" 
within the experimental chamber. The "wall" sepa- 
rating the two components is permeable to i (in the 
present case, i is hydrogen), but not to the other mixture 
components. Many materials are known to allow perme- 
ation of hydrogen. The best material for this purpose is 
palladium [23], a fact which has led to commercial ex- 
ploitation in hydrogen purifiers. Palladium itself has lit- 
tle mechanical strength, however, and in practical appli- 
cations, one must usually resort to the use of an alloy 
containing 25% silver. Palladium-silver, 75/25, is easily 
fabricated into useful forms such as rod, foil and tubing. 
Due to geometrical and support considerations, tubing 
is the most suitable form for use at high pressures. 

The Pd-Ag tubing was obtained from a commercial 
processor and had an outside diameter of 0.160 cm, and 
a wall thickness of 0.008 cm. The material is annealed to 
a dead soft condition in a tube furnace, and tested for 
pinholes using a helium leak detector. At present, the 
material is available only in 396 cm long sections, due to 
limitations of the drawing die. The thin wall of the 
tubing was necessary to achieve reasonable gas diffusion 
rates. Indeed, if the wall thickness were much greater, 
diffusion may not be observed at all [24]. It was neces- 
sary to increase the strength of this rather thin and weak 
tubing, since it must support pressure gradients of 7 MPa 
or more. For this purpose, 400 cm length coil springs 
were fabricated from 0.025 cm diameter stainless steel 
wire. The springs were wound to have a nominal re- 
laxed o.d. of 0.145 cm, and are easily inserted into the 
Pd-Ag tubes. The spring is then under a low tension 




Figure 3-Semipermeable membrane manifold. 

interference fit inside the Pd-Ag tube, and greatly in- 
creases the tube's "compressive load" strength. 

The Pd-Ag tubes each containing a coiled spring, 
were then formed into five layer coils, approximately 
2.5 cm long with a maximum outside coil diameter of 2.5 
cm. The coiling process was done on a custom made 
adjustable length mandrel held in a lathe collet. Extreme 
care was taken to keep all traces of dirt or oil from 
contaminating the tubing, since this could result in per- 
manent poisoning of membrane activity [23]. 

For each manifold (fig. 3), three coils were prepared 
and wound as described above. One end of each coil (the 
outside end) was capped off with a custom made 316 
stainless steel cap. The other end (the inside end) was 
fitted into a specially drawn "elbow" made from 0.208 
cm o.d. 321 stainless steel tubing, approximately 1.5 cm 
long. This elbow was, in turn, placed into a drilled hole 
in a central tube (0.328 cm o.d., 0.089 cm wall thickness). 
The central tube can be made from either 316 stainless 
steel or 270 nickel. It is the top of this central tube which 
is ultimately connected to the pressure measurement 
manifold using a compression fitting. The bottom of the 
tube is capped. 

The Pd-Ag coils, end caps, elbows, and central tube 
were assembled in a vibration eliminating fixture to al- 
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low brazing of all the joints. Due to the chemical sensi- 
tivity of Pd-Ag, the use of any kind of flux is impossible. 
Thus, vacuum brazing was employed. All joints were 
"dressed" with a single loop of 0.051 cm diameter pre- 
cious metal (82% gold, 18% nickel) brazing wire. The 
fixture-held manifold was placed in the vacuum furnace, 
which was then evacuated to 1.3 X 10 -5 Pa. The tem- 
perature was then raised to the melting point of the 
brazing wire, and the system was allowed to "soak" at 
this temperature for a few minutes. The furnace was 
then cooled slowly to prevent the development of 
cracks. The completed manifold was tested for pinholes 
using a helium leak detector, and for hydrogen diffusion 
(with pure hydrogen) on both sides of the membrane. 
Before the membrane is actually used, the Pd-Ag sur- 
face must be activated by repeated oxidation and reduc- 
tion at elevated temperature (300 °C has proven to be 
satisfactory). Oxidation is done in a pressurized (6 MPa) 
environment of oxygen or air; reduction is done in the 
same way using pure hydrogen [25]. A membrane pre- 
pared in this way will quickly reach equilibrium during 
an experiment. Typically, 45 minutes is adequate to al- 
low 95% of equilibration (pressure stabilization) to take 
place. 

3.3 Pressure Measurements 

Measurements of the pressures of the pure hydrogen 
(from inside the membrane manifold) and the mixture 
(outside the membrane) were made using a commercial 
dead weight pressure balance. The balance has an accu- 
racy of 0.015%, and a range of between 0.04 and 82.7 
MPa. The pressure of the test fluid (either hydrogen or 
mixture) is referenced to an inert gas line via a di- 
aphragm type differential pressure transducer. The inert 
gas communicates its pressure to another differential 
pressure transducer which is referenced to oil pressure. 
The oil pressure is generated by weights on a piston 
table. When both differential pressure transducers are at 
null condition, one can obtain the pressure of the test 
fluid by summing the weights pressurizing the piston. 
Corrections are made for local gravity (979,601.148 rail- 
ligals at Boulder station B), atmospheric pressure, tem- 
perature effect on piston area, elastic distortion of the 
cylinder, and oil head pressure. Corrections for fluid 
expansion into valve and transfer line volumes are also 
required. In actual practice, the differential pressure 
transducers are first brought to null with equal opposing 
pressures, to allow the setting of a zero reference point. 
Hysteresis effects in the transducer are minimized by 
"rocking" the diaphragm up and down several times 
before taking a measurement. The effect of temperature 
on the pressure transfer lines (between differential pres- 
sure transducers) is a small source of systematic error, so 
care was taken to maintain a uniform room temperature. 



4. Temperature Control and Measurement 

Since the fugacity measurements are made iso- 
thermally, adequate temperature control is necessary. 
Mixed liquid baths are usually the best way to thermo- 
stat an experiment of this type. This is impractical in the 
present work, however, since measurements at high 
temperatures are anticipated. For this reason, a commer- 
cially available vigorously mixed oven was modified to 
provide a thermostat. The major modifications included 
the installation of baffles (to promote more uniform mix- 
ing of the temperature transmitting fluid), installation of 
large aluminum thermal masses (to reduce temperature 
gradients) and the fabrication of a thermally insulated 
cover plate (from which the pressure vessel and other 
high temperature components are suspended inside the 
oven). 

The oven is heated by 6.5 kW resistive elements con- 
trolled by a proportional SCR circuit, allowing tem- 
perature control of the apparatus components to 
±0.05 "C. The controller is also equipped with a high 
temperature limiting cutoff for safety. For additional 
safety, the temperature transmitting fluid of the oven is 
maintained at a slight negative pressure, and a low flow 
rate of nitrogen is continually passed through the hot 
zone and is exhausted into a fume hood. 

Temperature gradients among key components inside 
the oven are monitored using opposed pairs of J-type 
(iron-constantan) thermocouples [26,27]. The observed 
gradients are then minimized using very low power 
"shimming" heaters located near or in major com- 
ponents. 

Temperature measurements are made using a com- 
mercial quartz crystal oscillator thermoprobe. The 
quartz probe is located in a thermowell in the pressure 
vessel. The immersion error of the probe in this config- 
uration was calculated to be on the order of ±0.003 °C. 
A single point calibration is performed on a regular 
basis, using either the triple point or freezing point of 
water. This single point measurement is for scaling only; 
the frequency versus temperature table for the quartz 
probe is stored on a ROM, and is accessed auto- 
matically. The probe has an absolute accuracy of 
±0.03 °C in the present region of interest. Errors due to 
long term instability amount to less than ±0.008 °C. 
Errors due to hysteresis are of a negligible level, since 
the measurements are made isothermally. Thus, ther- 
mometry precision on the order of ±0.05% (coefficient 
of variation) is easily obtainable. 

5. Analytical System 

It is necessary to determine the composition of the gas 
mixture (outside the semipermeable membrane) in order 
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to deduce fugacity coefficients, A custom designed de- 
velopmental gas chromatograph and sampling system 
was built to perform the necessary analyses. This was 
considered necessary because of some of the unique 
problems inherent in hazardous gas analysis [28,29]. 

The analytical system consists of a primary sampling 
system (which allows removal of a sample from the 
hostile pressure/temperature environment of the ves- 
sel), a secondary sampling system (which allows sample 
manipulation and injection) and the gas chromatograph, 
in which the separation and quantitation is actually 
done. The primary sampling system consists of a 10 port 
valve (rated for service at 250 °C and 48 MPa), a circu- 
lating pump, and necessary valving. This equipment is 
shown in place in figure 2. The sampling valve is actu- 
ated remotely using helium gas (at 0.83 MPa, feeding a 
pneumatic drive), thus permitting the operator to re- 
main outside the potentially hazardous pressure vessel 
zone. Helium is used as the actuation fluid to provide 
very fast valve switching. The evacuable sample loop of 
this valve has a volume of 0.5 ml. At a mixture pressure 
of 3.45 MPa, the withdrawal of a sample aliquot of this 
size causes a mixture pressure drop of 0.0008 MPa, 
which has a negligible effect upon the equilibrium mix- 
ture composition. Before the sampling valve is switched 
to remove this plug of sample, the mixture gas is circu- 



lated through the sample loop using a commercial check 
valve type pump, This minimizes the problem of con- 
centration gradients, which often plague expansions of 
gas mixtures from elevated pressures. Only a few 
minutes of circulation are needed, since the pump is able 
to move 4.0 ml of fluid per stroke, and a pumping rate of 
15 strokes per minute is employed. 

Switching of the primary sampling valve results in the 
transfer of the plug of sample to the secondary sampling 
manifold, which is shown schematically in figure 4. The 
sample is held between valves A and B inside a 316L 
stainless steel bellows. The pressure of the sample may 
be varied by compressing or expanding the bellows. 
Valve B allows the sample to be introduced into a 10 
port sampling valve, S, equipped with a 50 julI sample 
loop. Valve D allows the loop to be evacuated before 
being filled with sample and between analyses. It has 
been found in this work that loop evacuation before 
filling leads to a factor of two increase in precision of the 
measured chromatographic peak area counts. The vent 
valve, C, allows the sample pressure to be dropped to 
ambient before injection. The ambient pressure is mea- 
sured using a fixed cistern mercury barometer, with the 
appropriate corrections being applied to the readings 
[30]. 

The sampling valve, bellows, and valves A through D 
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are mounted within a massive (1050 cm 3 ) aluminum 
block which is heated (in the present work, to 
125±0.3 °C). The aluminum block reduces temperature 
gradients to a negligible level. Temperature control is 
provided by a proportional controller which responds 
to a thermistor sensor. Thus, by controlling the tem- 
perature and pressure of the fixed volume sample loop, 
one is assured of injecting the same quantity of sample 
for each analysis. 

Upon switching the sampling valve of the secondary 
manifold, the plug of sample is swept from the loop into 
the chromatograph. This chromatograph has been spe- 
cially constructed to provide precise control of column 
temperature and pressure. The column is thermostated 
in a vigorously stirred oil bath which is maintained to 
within ±0.02 °C. All analyses are done isothermally to 
minimize carrier flow rate disruptions. Column pressure 
is controlled using a fine adjustment pressure regulator 
(a commercially available unit used for pressure gauge 
calibration) followed by a fine metering valve. The 
chromatograph is very similar to units used for precise 
thermodynamic measurements [31-33]. Mass flow con- 
trollers are not used since the column temperature is not 
programmed. The carrier gas flow rate, measured at the 
detector exit, is maintained at 40.00+0.05 ml per 
minute. Since both hydrogen and the mixture gas (in the 
present case, methane and propane were studied) re- 
quired detection and quantitation, argon was chosen as 
the carrier gas [34]. 

The separation of mixtures of hydrogen-methane and 
hydrogen-propane is experimentally trivial and will be 
described later. A micro cell thermal conductivity de- 
tector (TCD) at the column exit, connected to a com- 
mercial electronic integrator, provided detection and 
quantitation. Peak area counts as logged by the integra- 
tor are corrected for atmospheric pressure (equivalent 
to sample loop pressure, since injections were done at 
ambient pressure) and detector aging. The pressure cor- 
rection is required since the loop pressure affects the 
quantity of fluid injected, and dividing by the atmo- 
spheric pressure normalizes the response. The cor- 
rection for detector aging is needed since the gold- 
sheathed tungsten filaments of the TCD run hotter in an 
argon carrier stream than in the more common carrier 
gases such as helium and hydrogen. Thus, the filaments 
undergo accelerated oxidation due to trace amounts of 
oxygen that are not trapped. This correction is applied 
by periodically injecting a standard mixture to track 
signal loss. The corrected peak areas are then converted 
to mole fractions either by using a calibration equation 
or the standard bracketing technique [35,36]. After leav- 
ing the detector exit, sample and carrier are expelled 
through a transfer line to a dedicated outside blower, to 
provide operator safety. 



6. Testing and Evaluation 

Since the analytical measurements are the greatest 
single source of error in this experiment, a thorough 
testing of the sampling and chromatographic systems 
was needed. Testing of the nonanalytical components of 
the apparatus was described in the preceding section. In 
this section, the extensive tests performed on the anal- 
ytical system will be described. 

Due to the complexity of the two sampling systems, 
with many joints, fittings and valves, correlation be- 
tween repeated measurements is very possible. For ex- 
ample, if sampling valves are not sufficiently evacuated 
or purged, residual sample from previous analyses 
would cause a systematic, additive error. This would be 
reflected in a definite trend of component areas with 
analysis repetition or time. 

To check for this type of sample correlation, an ex- 
tended study was done on a standard mixture of 74.8% 
hydrogen, 25.2% methane. The mixture was prepared as 
described (for standard calibration mixtures) in the next 
section. Thirty separate analyses were done on the mix- 
ture over a period of two days. The measurements were 
made in approximately equal intervals of time. Com- 
ponent areas, total areas, and component chro- 
matographic mole fractions were calculated. Chro- 
matographic mole fractions are those determined from 
the component area and total area, uncorrected for rela- 
tive detector response. It was then necessary to look for 
trends in these measured quantities as a function of sam- 
ple run or time. 

Plots of each of the above quantities versus run (or 
time) revealed no discernible structure. In addition, ta- 
bles of correlation were calculated, which included 
Pearson, Spearman, and Kendall correlation coeffi- 
cients. No unexpected correlations were observed. Re- 
lated variables, such as component area and total area, 
showed high correlation as expected [37,38]. Linear re- 
gressions were performed in which the measured chro- 
matographic quantities were fit against the repetition 
number. In each case, little correlation was found. The 
estimates of the first order parameter were poorly deter- 
mined (as deduced from the standard errors) and of 
questionable significance (as deduced from the t values), 
while the intercepts were well determined and highly 
significant. This suggests that the proper representation 
for repeat analysis data is a mean with some measure of 
dispersion. Since the data were taken at approximately 
equal intervals of time, time,series statistics are nomi- 
nally informative. Thus, first order auto-correlation and 
Durbin- Watson d statistics were calculated, and indi- 
cated minimal correlation among errors [39]. Since an 
average appears to represent repetitive analysis ade- 
quately, we can make probability plots to assess the 
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nature of the distribution of data about the mean. These 
tests showed that the actual distribution of measurement 
values deviates negligibly from the normal distribution, 
and inferences based on the assumption of normality are 
valid [40]. 

It may be concluded from the results of foregoing 
tests that each analysis is an independent measurement 
of concentration, unaffected by previous measurements, 
and that an average of repetition measurements will pro- 
vide an unbiased estimate of concentration, and the stan- 
dard deviation will provide a measure of the dispersion. 



conditions as the experimental mixtures (same condi- 
tions of carrier flow rate, injector, column and detector 
temperature, and TCD bridge current). The standard 
mixtures were prepared gravimetrically using a com- 
mercial 25 kg capacity two pan balance. The sensitivity 
of this balance is ±0.0025 g, with an accuracy de- 
pendent on calibration with a standard 1 g class S weight 
(1.0000±0.0005 g) [41]. This uncertainty will propagate 
an error of approximately 0.4% into a typical gravimet- 
rical mole fraction (a mixture of 0.5 mole fraction hydro- 
gen prepared by three separate weight measurements). 



7. Quantitation 

Since the sensitivity of the thermal conductivity de- 
tector (TCD) is sample material dependent, it is not 
possible to take the integrator response as a concen- 
tration reading. It should be noted that absolute de- 
tectors requiring no calibration are available (for exam- 
ple, a gas density balance), and for these instruments 
direct integrator response is a valid method of quan- 
titation. For relative response detectors such as the 
TCD, one must use either an internal or external stan- 
dard, or normalization [35,36]. For this type of work, 
the external standard method is preferable. 

To use the external standard method, mixtures of 
known composition of hydrogen/methane and 
hydrogen/propane were prepared. These standard gas 
mixtures were then chromatographed under the same 



8. Hydrogen-Methane Calibration 
and Analysis 

Nine standard mixtures of hydrogen-methane were 
prepared for the calibration, along with a sample of pure 
hydrogen. Pure methane was not used since it represents 
an experimentally unaccessible point in the mole frac- 
tion range. Separation was easily done using a packed 
column (1.5 m long, 0.32 cm o.d.) of 150-200 mesh 5 A 
molecular sieve. The column was maintained at 60 °C, 
since methane tends to adsorb on the zeolite surface of 
the molecular sieve at lower temperatures. The elevated 
temperature also provided for increased speed of anal- 
ysis without sacrificing baseline resolution of the peaks. 
A plot of the corrected area counts (detector response) 
versus mole fraction is provided in figure 5. The slight 
but noticeable curvature of this plot is not unusual for 
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hydrogen-containing binary mixtures. 

The data were fitted to a second degree polynomial, 
with the result that the poorly determined constant term 
proved to be statistically indistinguishable from zero. 
This is internally consistent, since the curve must log- 
ically pass through the origin, A second degree poly- 
nomial constrained to pass through the origin (i.e., the 
constant term is set equal to zero) was then found to 
represent the data to within experimental error. The 
results of the analysis are shown in table 1. The resulting 
parameters a and b were highly significant (as judged by 
the /-test and the ^-test), although a is far better deter- 
mined than is b. The percent error in the fit is 0.55, 
which is comparable to that obtained using other chro- 
matographic quantitation methods. The root mean 
square error, or overall regression standard deviation, 
was ±0.0026 in mole fraction, resulting in an error of 
0.52% in mole fraction for an equimolar mixture. A plot 
of residuals versus mole fraction showed no discernible 
structure, and a regression found no meaningful cor- 
relation. Thus, the model set forth in table 1 accounts for 
all data structure except random error. 

Table 1. Regression results for hydrogen-methane calibration 



aA -\-bA 2 =n 

where A — corrected chromatographic area response 

n = mole Fraction of hydrogen 

regression results: 

a = 8 .5764 x 10" ' 0.56 percent CV 
b= -2.5940 x 10-' 4 38.70 percent CV 
r 2 =0.99998 

root mean square error, n =0.0026 
coefficient of variation =0.5506 



9. Hydrogen Propane Calibration 
and Analysis 

Standard calibration mixtures were prepared for 
hydrogen-propane in the same way as for hydrogen- 
methane. For this system, only four standard mixtures 
were prepared, and only in the higher hydrogen mole 
fraction range. Calculations were done using the stan- 
dard "bracketing" method [29] rather than by deter- 
mining a calibration equation. The reason for this type 
of calibration (rather than the more preferable type done 
for hydrogen-methane) is the very high cost of research 
grade propane. It is, however, not unreasonable to ex- 
pect mole fraction accuracy of better than 1.0% for the 
hydrogen-propane system. 



This mixture is easily separated using a packed col- 
umn (1.5 m long, 0.32 cm o.d.) containing porapak-Q 2 . 
The column was maintained at 80 °C, since the porous 
polymer packing was observed to entrain hydrogen at 
lower temperatures.. It should be noted that elevated 
temperatures cause the polymer beads to swell, resulting 
in an increase in the retention time of hydrogen and a 
decrease of that of propane. Thus, the temperature at 
which the two peaks merge is lowered by these additive 
effects. 



10. Results for Hydrogen/Propane and 
Hydrogen/Methane 

Plots of the fugacity coefficient of hydrogen in mix- 
ture with propane, <f>H 2 , as deduced from eqs (2), (3), and 
(4), are shown in figures 6 and 7. The data shown in 
figure 6 were taken at 80 °C (353 K), and those in figure 
7 were taken at 130 "C (403 K). The nominal mixture 
(total) pressure in each case was 3.45 MPa. The actual 
data are listed in tables 2 and 3. The error bars represent- 
the uncertainty of a typical <$>h 2 — *h 2 pair, and will be 
discussed in more detail later. The data represented by 
the triangles on figure 6 are the results of an earlier study 
[17]. Although the earlier data were taken at a slightly 
different temperature (345 K instead of the 353 K in the 
present work), the agreement is striking. 

In both plots, the change in 4>h 2 is most pronounced at 
lower values of x Hl . This behavior is consistent for a low 
molecular weight gas (at a relatively high reduced tem- 
perature) in a binary mixture with a heavier gas. It 
should also be noted that the data at 1 30 °C show a more 
gentle increase in slope than the data at 80 °C. This is 
consistent with the expectation that <J>h 2 be closer to 
unity at higher temperatures, for a given composition. 
The slight upturn at the high hydrogen mole fraction 
range is probably due to the upturn in molar volume in 
this region. This behavior of <J>h 2 is also suggested in an 
earlier study done under approximately the same condi- 
tions [15]. A study comparing the measured values of 
<J>h 2 with predictions of several equations of state is on- 
going and will be reported later [42]. 

A plot of 4>% 2 for hydrogen/methane at 80 °C (353 K) 
and 3.45 MPa is provided in Figure 8. The actual data are 
listed in table 4. The total change of <J>h 2 over the hydro- 
gen mole fraction range studied is 13%, compared to 



2 Certain commercial equipment, instruments, or materials are iden- 
tified in this paper in order to adequately specify the experimental 
procedure. Such identification does not imply recommendation or 
endorsement by the National Bureau of Standards, nor does it imply 
that the materials or equipment identified are necessarily the best 
available for the purpose. 
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Figure 6-Fugacity coefficient of 
hydrogen, (f>fi 2 , versus mole frac- 
tion of hydrogen, in propane at 
80 °C (353 K). A typical error 
bar is shown at the lower left. 
o — this work 

^ — previous data of Antezana 
[17], taken at 345 K. 
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Figure 7-Fugacity coefficient of d> r 
hydrogen, tj>tf 2 , versus mole frac- 
tion of hydrogen, in propane at 
130 "C (403 K). A typical error 
bar is shown at upper right. 
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23% in the hydrogen-propane mixture at 353 K. This is 
expected because of the higher reduced temperature of 
methane. The value of <f>Hj changes most sharply in the 
region of low hydrogen mole fraction, as expected [10]. 
Data on this system were not taken above 353 K since 
the change in </>% would probably not be detected. As 
with the hydrogen-propane system, work on predictions 
and comparison with equation of state results is in 
progress and will be reported on later. 

The error bars shown in figures 6-8 represent random 
propagated errors in 4>h 2 and x Hl . From eq (2), only the 



first term on the righthand side was considered, since 
the contribution of the second order term is small. Er- 
rors in temperature, hydrogen partial pressure, and the 
second virial coefficient were then considered. Strictly, 
the error in B might properly be considered systematic, 
however enough data are available on hydrogen to as- 
cribe a measure of dispersion about B [43-45]. The error 
propagation for eqs (3) and (4) was carried out in the 
usual manner [46]. The recognized sources of systematic 
error, along with countermeasures, have been described 
in the experimental section. 
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Table 2. Hydrogen/propane at 80 °C (353 K) 



Table 4. Hydrogen/methane at 80 °C (353 K) 



Mole fraction, Hj 


Partial pressure, H2 
(MPa) 


Fugacity coefficient, <f>ff 2 


Mole fraction, H2 


Partial pressure, H 3 
(MPa) 


Fugacity coefficient, 4>Si 


0.848 


3.01 


1.040 


0.849 


3.17 


1.121 


0.822 


2.91 


1.033 


0.790 


2.92 


1.115 


0.730 


2.73 


1.028 


0.764 


2.80 


1.097 


0.659 


2.33 


1.036 


0.652 


2.40 


1.115 


0.594 


2.15 


1.058 


0.536 


1.99 


1.125 


0.526 


1.96 


1.080 


0.459 


1.75 


1.134 


0.445 


1.69 


1.106 


0.307 


1.19 


1.166 


0.344 


1.39 


1.170 


0.216 


0.88 


1.223 


0.309 
0.280 


1.33 
1.25 


1 230 








1.283 









Table 3. Hydrogen/propane at 130 "C (403 K) 



Mole fraction, Hj 


Partial pressure, H2 
(MPa) 


Fugacity coefficient, <f>K 2 


0.904 


3.29 


1.058 


0.897 


3.23 


1.054 


0.835 


2.96 


1.038 


0.783 


2.81 


1.047 


0.738 


2.66 


1.047 


0.618 


2.28 


1.077 


0.545 


2.06 


1.096 


0.472 


1.82 


1.116 


0.265 


1.12 


1.220 



The author would like to thank Drs. Gerald C. Straty 
and Howard J. M. Hanley for many helpful discussions, 
and Gretchen L. Hume for her assistance in per- 
forming some of the measurements. 
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